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IMPROVED STRESS-IMTEHSITY FACTORS FOH SEMI-ELLIPTICAL 
SURFACE CRACKS IN FINITE-THICKNESS PLATES 


I. S. Raju and J. C. Newian9 Jr. 

NASA-Langley Research Center 
Hampton, Virginia 

SUMMARY 

Surface cracks are among the nore common Vlavs in aircraft and 
pressure vessel components. Accurate stress analyses of surface-cracked 
components are needed for reliable prediction of their crack-grovth 
rates and fracture strengths. Several calculations cf stress-intensity 
factors for semi-elliptical surface cracks subjected to tension have 
appeared in the literatiire. However, some of these solutions are in 
disagreement by 50 to 100 percent. 

In this paper stress- intensity factors for shallow and deep semi- 
ellipticeil surface cracks in plates subjected to tension are presented. 
To verify the accuracy of the three-dimensional finite-element models 
employed, convergence was studied by varying the number of degrees of 
freedom in the models from 1500 to 6900. The 69OO degrees of freedom 
used here were more than twice the number used in previoivsly reported 


solutions. 
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INTRODUCTION 

t'urface cracics are among the more cosmion flaws in aircraft and 
pressiire vessel components. Accurate stress analyses of these surface- 
cracked ccMsponents are needed for reliable prediction of their crack- 
grovtl: rates and fracture strengths. Exact solutions to these difficult 
problems are not available; therefore, approximate methods must be used. 
For a semi-elliptical surface crack in a finite-thickness plate subjected 
to tension (Figure l). Browning and Staiith [l], Kobayashi [2], and Smith 
and Sorensen [3] used the alternating method and Kathiresan [4] used 
the finite-element method to obtain the stress- intensity factor varia- 
tions along the crack front for various crack shapes. For a de^ semi- 
elliptical surface crack (with a/t * 0.8 and a/c = 0.2) subjected to 
tension, the stress-intensity factors obtained by Smith and Soroisen [3]» 
Kobayashi [2], and Kathiressin [U] disagreed by 50 to 100 j>ercent. The 
reasons for these discrepancies are not well understood. 

This paper presents stress-intensity factors for shallow and deep 
semi-elliptical surface cracks in plates subjected to uniform tension. 

To test the validity of the present analysis, two crack configurations, 
both embedded in a large body subjected to uniform tension, were analyzed: 
(l) a circular (penny- shaped) crack eind (2) an elliptical crack. These 
results are compared with exact solutions from the literature [5l» To 
verify the accuracy of the solutions for surface cracks in finite- 
thickness plates, convergence was studied by varying, from 1500 to 69OO, 
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the number of degrees of freedom in the finite-elonent models. The 
6900 degrees of freedom used here were more than twice the largest 
number used previously [b]. These models were cooqposed of singular 
elements around the crack front and isoparametric (linear strain) 
elements elsewhere. Mode I elastic stress- intensity factors were 
calculat i by using a nodal force method. 

The computations reported herein were conducted on a unique 
computer called the STAR-100 at the NASA Langley Research Center. 
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THR]2E-DIMENSI0IiHL ANALYSIS 

Several a^ tempts [6-8] have been made to develop special three- 
dimensional finite-elements that account for the stress and strain 
singularities caused by a crack. These elements have had assumed 
displacanent cr stress distributions that simulate the square-root 
singularity of the stresses and strains at the craek front. Such 
three-dimensioned singularity elements [ 9 ] were also used herein for 
the analysis of finite-thickness plates containing embedded elliptical 
or sani-elliptical surface cracks (see» for example. Figure l). 

Finite-Element Idealization 

Two types of elements (isoparametric and singxdar) [ 9 ] were used 
in combination to model elastic bodies with embedded elliptical cracks 
or s«tti-elliptical surface cracks. Figure 2 shows a typical finite- 
element model for an embedded circular crack in a large body 
(h/a = b/a = 5 )« This model idealizes one eighth of the body. Various 
numbers of wedges were used to form the desired configuration. Fig- 
ure 2(a) shows a typical model with eight wedges. Each wedge is 
composed of elements that are identical in pattern to that shown in the 
(J) = constant plane. The arrangement of the elements aroiind the crack 
front is shovn in Figure 2(b). The isoparametric (linear strain) 
elements (denoted as I) vere used everywhere except near the crack 
front. Around the crack front each wedge contained eight ** singularity’* 
elements (S) in the shape of pentahedrons. The "singularity" elenents 
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had square-root terms in their assumed displac«nent distribution and, 
therefore, produced a singular stress field at the crack front. Details 
of the formulation of these types of elements are given in Reference 9- 
'Phe finite-element model for the embedded elliptical or semi- 
elliptical surface crack was obtained from the finite-el«nent model for 
the circular crack by using an elliptic transformation. This transfor- 
mation was needed because the stress- intensity factors must be evaluated 
from either crack-opening displacements or nodad forces along the normals 
to the crack front. If (x,y,z) are the Cartesian coordinates of a node 
in the circular-crack model and (x*,y*,z*) are the coordinates of that 
same node in the elliptical-crack model, then the transformation is 
given by 



for X and z not at the origin. Figure 3 shows how circular arcs and 
radial lines in the x,z plane of the circiilar-crack model are trauis- 
formed into ellipses and hyperbolas, respectively, in the x*,z* plane 
of the elliptical-crack model using equations (l). Because equations (l) 
are not vedid at the origin, a circle of very small radius, a/lOOO, was 
used near the origin in the circular crack model. The small circle 
maps onto an extremely narrow qfllipse in the x’,z* plane. The use of 
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the small circle avoids ill-shaped elements near the origin in the 
elliptical-crack model. Figure k shows a typical finite-element model 
of a finite plate containing ar elliptical crack. The transformation 
reduced the b/c ratio; therefore, in order to maintain b/c ^ 1*, 
additional elements were added along the x*-axis to eliminate the 
influence of plate width. 

Stress- Intensity Factor 

The stress-intensity factor is a measure of the magnitude of the 
stresses near the crack front. Under general loading, the stress- 
intensity factor depends on three basic modes of deformation (tension 
and in- and out-of-plane shear). But here only tension loading was 
considered and, therefore, only Mode I deformations occurred. The 
Mode I stress-intensity factor, K, at any point along an elliptical 
or semi-ellipticeLL crack in a finite plate is taken to be 

K = S F(f. i 4.) (2) 

where S is the applied stress, a is the crack depth, Q is the shape 
factor for an ellipse and is given by the square of the complete 
elliptic integral of the second kind [5]. The boundary-correction 
factor, F (a/t, a/c, 4>), is a function of crack depth, crack length, 
plate thickness and the parametric angle of the ellipse. The present 
paper gives values for F as a function of a/t and (\> for a/c = 0.2 and 
1.0. The a/t values ranged from 0.2 to 0.8. 
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The stress-intensity factors from the finite-element models of 
Quhedded elliptical and semi-elliptical surface cracks were obtained 
by using a nodal force method, details of which are given in Reference 9* 
In this method, the nodal forces normal to the crack plane and ahead of 
the crack front are used to evaluate the stress- intensity factor. In 
contrast to the crack-opening displacement method [6], this method 
requires no prior assumption of either plane stress or plane strain. 

For a surface crack in a finite plate, the state of stress varies from 
plane strain in the interior of the plate to plane stress at the 
surface. Thus, the crack-opening displacement method coxxld yield 
erroneous stress- intensity factors. 

RESULTS AND DISCUSSION 

In the following sections a circular (penny- shaped) crack and an 
elliptical crack completely embedded in a large body su'iected to 
uniform tension were analyzed using the finite-element method. The 
calculated stress-intensity factors for these crack configurations are 
compared with the exact solutions to verify the validity of the present 
finite-element method. 

For a semi-ci: ,ular and semi-elliptical surface crack in a finite- 
thickness plate, convergence of the stress- intensity factors was studied 
while the number of degrees of freedom in the finite-element models 
ranged from about 1500 to 6900. The stress-intensity factor variations 
along the crack front for semi-circular (a/c = 1) and semi-elliptical 
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(a/c = 0.2} surface cracks were obtained as functions of a/t vlth 
h/c = b/c ^ 4 . Whenever possible, these stress-intensity factors are 
compared with results from the literature. 

Exact Solutions 

In this section a comparison is made between the stress- intensity 
factors calculated from the finite-elenent analysis and from the exact 
solution [ 5 ] for an embedded circvdar (penny-shai>ed) crack (a/c = l) 
and an embedded elliptical crack (a/c = 0.2) in an infinite body. In 
the finite-element model, h and b were taken to be large enoiigh that 
the free boundary would have a negligible effect on stress intensity. 

The boundary correction on stress intensity for a circulsu* crack in a 
cylinder of radius b with b/a = 5 is about one percent [lO]. Therefore, 
to simulate a large body the finite-element model was assigned the 
dimensions h/a = b/a = 5, along with 3078 degrees of freedom. The 
calculated stress-intensity factors along the crack front from the 
circulau: crack model were about O.lt percent below the exact soli 'ion. 

The embedded elliptical crack (a/c = 0.2) model, with h/c = h/c - k 
and t/a = 5, had 33^8 degrees of freedom. The influence of finite 
boundaries on stress intensity was estimated to be about one percent. 

The finite-element analysis gave stress-intensity factors along the 
crack front generally within one percent of the exact solution, except 
in the region of sharpest curvature of the ellipse, where the calculated 
values were about 3 percent higher than the exact solution. Further 
refinement in the mesh size in this area gave more accurate stress- 
intensity factors. 
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Because the present method yielded stress- intensity factors for 
completely embedded circular and elliptical cracks within 0.4 to 
3 percent of the exact solutions, the method was considered suitable 
for analyses of more complex configurations, provided that enough 
degrees of freedom were used to obtain good convergence. 

Approximate Solutions 

To verify that the finite-element meshes used for the circular 
and elliptical crack models were sufficient to analyze cracked plates 
with free boundaries, through-the-thickness cracks with crack length- 
to-width ratios ranging from 0.2 to 0.8 were analyzed under plane 
strain assumptions. The meshes used here were exactly the same as 
those which occur on the (j) = tt/ 2 plane of the circular and elliptical 
crack models. These meshes were then used to model the cer.ter-crack 
tension specimen. For crack length-to-width ratios (c/b) ranging from 
0.2 to 0.6, the finite-element results were within 1.3 percent of the 
approximate solutions given in Reference 10. For c/b = 0.8, the finite- 
element result was 2 percent below the solution given in Reference 10. 
These results indicate that the mesh pattern along any plane has 
enough degrees of freedom to accotint for the influence of free 
boundaries on the stress- intensity factor. 

Convergence 

In the previous section the mesh pattern along any (J) = constant 
plane was found to be sufficient to account for the influence of free 
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boundaries even for very deep cracks. In this section, the mesh 
pattern in the angular direction, (|), is studied. Figures 5 and 6 show 
the results of a convergence study on the stress-intensity factors for 
a semi-circular surface crack and a s ^mi-elliptical surface crack, 
respectively, in a finite-thickness plate. The irger numbers of degrees 
of freedom are associated with smaller wedges in the (t-direction and 
a more accurate representation of the cruck shape. 

For the semi-circular crack (Figure 5)» h/a = b/a - 5 and 
a/t = 0.8. This configuration was chosen because the close proximity 
of the back siarface to the crack front was expected to cause difficulty 
in achieving convergence. The munber of degrees of freedom ranged from 
1500 to 6195 • The two finest models (U 3 IT and 6195 degrees of freedom) 
gave stress-int entity factors within about one percent of each other. 
Therefore, the moC il with U317 degrees of freriom was used subsequently 
to obtain stress- intensity factors as a fiinci ion of o/x. 

Figure 6 shows the convergence study for a semi-elliptical 
surface crack ^ ,c = 0.2). The h/c and b/c ratios were equal to four 
and, again, a/t was chosen as 0.8. The number of degrees of freedom 
ranged from xo92 to 6867 . The model with ^797 degrees of freedom gave 
results within about one percent of those from the finest model. 
Therefore, the model with ^797 degrees of freedom was used subsequently 
to obtain stress-intensity factors and crack-opening displacements as 
a function of a/t. 
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SeKi-Circular Surface Crack in a Finite-Thickness Plate 
Figure 7 shows and Table 1 presents the stress-intensity factors 
for a semi-circular surface crack in a finite-thickness plate as a 
function of the parametric angle, and the crcu:k depth- to-plate 
thickness ratio, a/t. Bear the intersection of the crack and the free 
surface, the stress-intensity factor increases more rapidly with a/t than 
at the deepest point ($ = v/2). For each value of a/t, the .'tress-intensity 
factor calculated from the model with 1»317 de yees of freedom is largest 
at the free surfeu:e = 0). 

For a semi-circiJ.ar surface crack with a/t = 0.55, the stress- 
intensity factors calculated by Browning and Smith flj, who used the 
alternating method, were about 1 to 3 percent below the present results 
for various values of <p. 

Semi-Elliptical Surface Crack in u Finite-Thickness Plate 
Figure 8 shows and Table I prosents the stress-intensity factors 
for a semi -elliptical s'orface crack (a/c = 0.2) in a finite-thickness 
plate as a function of the parametric angle, (^, and the crack iepth-to- 
plate thickness ratio, a/t. For each value of a/t, the maximum stress- 
intensity factor occurs at the deepest point = n/2). Also, the maxim»un 
stress- intensity factor is larger for larger values of a/t. 

Figures 9 and 10 show stress-intensity factors obtained by several 
investigators for a semi-elliptical surface crack in a finite-thickness 
plate. Figure 9 shows the results fer a surface crack with a/c = 0.2 
and a/t = 0.8. Smith and Sorensen [3] and Kobayashi 12] used the 
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alternating method and Kathiresan [fe] used the finite-element method 
to obtain stress-intensity factor variations along the crack front. 
These three solutions disagree by as much as 50 to 100 percent. The 
reasons for these discrepancies are not veil understood. The present 
results, shown as solid symbols, are considerably higher than the 
previous solutions [2-k]. The results from Staith and Sorensen [3] are 
generally closer to the present results, though 10 to 25 percent lover. 

Figure 10 shows a comparison of the mazimum stress-intoisity 
factors obtained by several investigators for a seni-elliptical surface 
crack as a function of a/t. The maximim stress-intensity factors 
occurred at 4> = v/2. The present results are shown as solid symbols. 
The open sjrmbols show the results from anith and Sorensoi [3l* 

Kbbayashi [2,U], and Rice and Levy [l2]. The results fz'om Rice and 
Levy, obtained from a line-spring model, are about 3.5 percent below 
the present results over an e/t range from 0.2 to 0.6. For a/t > 0.6, 
the Rice-Levy solution shows a reduct i<» in stress intensity. The 
dash-dot curve shows the results of an approximate equation proposed 
by Newman [13] for a vide range of a/c and a/t ratios. Newman's 
equation is. within ±5 percent of the present results over an a/t r. nge 
from 0.2 to 0.8. Newman's equation gives a good engineering estima.e 
for the maximum stresii-intensity factor. 
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CONCLUDING REMARKS 

A three-dijaensional finite-element elastic stress inalysis vas 
used to calculate stress- intensity factor variations along the crack 
front for completely ^bedded elliptical cracks in l€urge bodies and 
for semi-elliptical surface cracks (crack depth-to-crack length ratios 
vere 0.2 and l.O) in finite-thickness plates. Three-dimensional 
singularity elu&ents vere used at the crack front. A nodal force 
method irtiich requires no prior assumption of either plane stress or 
plane strain vas used to evaluate the stress-intensity factors along 
the crack front. 

Completely onhedded circular emd elliptical cracks vere analyzed 
to verify the accuracy of ^he finite-element analysis. The stress- 
intensity factors for these cracks vere generally about O.U to 1 percent 
belov the exact solutions. Hovever, for the elliptical crack the 
calculated stress-intensity factors in the region of sharpest curvature 
of the ellipse vere about 3 percent higher than the exMt solirtion. 

The numbers of degrees of freedom in the embedded crack models vere 
about 3000. A convergence study on stress-intensity factors for semi- 
ellipticcd. surface ciacks in finite-thickness plates shoved that 
convergence vas achieved for both the sani-circular and the semi- 
elliptical surface crack vith about U 5 OO degrees of freedom. 

For the semi-circular surface crack the maxinam stress- intensity 
factor occurred at the intersection of the crack vith the free surface. 
On the other hand, for the semi-elliptical sxjrface crack (crack depth- 
to-crack half length ratio of 0.2), the maximum stress- intensity factor 
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occxirred at the deepest point. For both the semi-circular and semi- 
elliptical surface crack the stress-intensity factors vere larger for 
larger values of crack depth-to-plate thickness ratio. 

For the semi-circular surface crack* the stress-intensity factors 
calculated by Browning and Smith Ll] using the alternating method, 
agreed generally within about 3 percent with the present results. 
However, for semi-elliptical s\irface cracks (crack depvh-to-plate 
thickness ratio of 0.2) Snith and Sorensen [3l using the alternating 
method gave stress-intensity factors in considerable disagre^nest 
(10 to 23 percent) with the present results. For semi-elliptical 
surface cracks the results frcm Rice and Levy [12] for crack d^th- 
to-plate thickness ratios less than or equal to 0.6 and an approximate 
equation proposed by Newman [13] vere in good agreement with the 
present results. 

The stress- intensity factors obtained herein shoxild be useful in 
correlating fatigue crack-growth rates as well as fracture toughness 
calctilat ions for the surface-crack configuration.' considered. 
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TABLE I: Boundary Correction Factors, F, for Semi-Elliptical 

Surface Cracks Subjected to tension 
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(V = 0 . 3 ; F = K/S /ira/Q) 
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(a) Circular crack. (b) Elliptic crack. 


Fig. 3.- Circle to ellipse transformation. 
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Fig. 7.- Distribution of stress-intensity factors along crack front 
for a semi-circular surface crack (Q = rr^/A). 







Fig. 10.- Comparison of the maximum stress-intensity factor for a 
semi-elliptical surface crack as a function of a/t 
(^= 7T/2; Q = 1.104; a/c = 0.2). 


